Morphologies and flow patterns in quenching of lameUar systems with shear 
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We study the behavior of a fluid quenched from the disordered into the lamellar phase under the 
action of a shear flow. The dynamics of the system is described by Navier-Stokes and convection- 
diffusion equations with pressure tensor and chemical potential derived by the Brazovskii free-energy. 
Our simulations are based on a mixed numerical method with Lattice Boltzmann equation and finite 
difference scheme for Navier-Stokes and order parameter equations, respectively. We focus on cases 
where banded flows are observed with two different slopes for the component of velocity in the 
direction of the applied flow. Close to the walls the system reaches a lamellar order with very few 
defects and the slope of the horizontal velocity is higher than the imposed shear rate. In the middle 
of the system the local shear rate is lower than the imposed one and the system looks as a mixture 
of tilted lamellae, droplets and small elongated domains. We refer to this as to a region with a Shear 
Induced Structures (SIS) configuration. The local behavior of the stress shows that the system with 
the coexisting lamellar and SIS regions is in mechanical equilibrium. This phenomenon occurs, at 
fixed viscosity, for shear rates under a certain threshold; when the imposed shear rate is sufficiently 
large, lamellar order develops in the whole system. Effects of different viscosities have been also 
considered: The SIS region is observed only at low enough viscosity. We compare the above scenario 
with the usual one of shear banding. In particular, we do not find evidence for a plateau of the 
stress at varying imposed shear rates in the region with banded flow. We interpret our results as 
due to a tendency of the lamellar system to oppose to the presence of the applied flow. 

PACS numbers: 64.75.+g; 05.70.Ln; 47.50.-|-d; 82.35. Jk 



I. INTRODUCTION 

Complex fluids such as polymer solutions, liquid crys- 
tals, or surfactant systems are characterized by the pres- 
ence of organized structures at mesoscopic scales between 
macroscopic and solvent molecular lengths 0. Under 
the action of external forcing, the coupling between the 
mesoscopic structures and the local velocity field makes 
the flow properties of complex fluids different from those 
of simple liquids ^211. For example, when a shear flow is 
applied to a simple fluid, after a short initial transient, 
the linear relation between the stress S and the shear 
rate 7 is verified, 5* = 777 with ij being the fluid viscos- 
ity. In complex fluid an effective viscosity can be deflned 
by the same relation but its value is not constant, de- 
pending on the strength of the applied flow. In systems 
with interfaces, for examples, the effective viscosity gen- 
erally decreases (shear thinning) when the shear rate is 
increased due to the alignment of interfaces with the flow. 

In some cases the flow can induce new organization in 
the fluid '35. In lamellar phases, for instance, onion and 
hexagonal phases have been observed not existing at rest 
UilEP' Such shear induced structures (SIS) can coexist 
in a range of applied shear rates with the structures un- 
modified by the flow. Usually the SIS and the unmodified 
phase have different viscosities so that they flow with dif- 
ferent profiles. This phenomenon, called shear banding. 



has been observed in many complex fluids and also in 
lamellar systems 0; • Its microscopic origin is not yet 
completely understood and different explanations have 
been proposed depending on the system Q • 

The traditional theoretical description of shear band- 
ing and other rheological behaviors in complex fluids is 
based on the assumption of a local relation between the 
stress and the shear rate However, in this way, 

the role of the structures present in the fluid is not evident 
and they cannot be directly related to the flow pattern. 
While a full description of the system with all its molec- 
ular dynamical variables is not possible, a description at 
mesoscopic level based on an order parameter evolution 
equation can enlighten many issues concerning the ki- 
netics of shear banding, the morphology of the different 
structures in the fluid and the evolution of the flow field 
[Tl. IT^ . Since the formation of SIS often occurs during 
transient regimes and since, also under stationary con- 
ditions, the flow pattern is not known a priori, one un- 
derstands that a description based on both Navier-Stokes 
and order parameter equations is generally required. The 
hydrodynamical description is useful also for comparisons 
with experiments where not always the dynamical quan- 
tities are all easily accessible. 

The aim of this paper is twofold. First, we want to 
show the relevance of the full hydrodynamical description 
also in systems with applied flow and see how a speciflc 
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problem can be conveniently studied by our simulation 
methods. Then we hope to clarify some aspects of the 
kinetics of formation of lamellar phases in cases when 
banded flows are observed 14]. In our model the lamel- 
lar properties are encoded in a free-energy functional of 
an order parameter representing the relative concentra- 
tions of substances in the mixture. The order parameter 
follows a convection-diffusion equation; from its instan- 
taneous configurations the stress can be calculated and 
inserted in the Navier-Stokes equation without assuming 
a further stress-shear rate constitutive relation. 

The free-energy functional considered in this paper 
was originally introduced to study the effects of fluctua- 
tions in the disordered-lamellar phase transition [l5l | and 
later used to describe equilibrium properties of di-block 
copolymer systems Solutions of copolymers, con- 

sisting of A-polymers covalently bonded to B-polymers 
in pairs, can organize in striped phases where A-rich and 
B-regions are separated by a stack of lamellae ■ The 
order parameter in this case represents the relative con- 
centration of A and B substances. Our model is also rel- 
evant for other systems with lamellar order. We mention 
ternary mixtures where the surfactant form interfaces be- 
tween oil and water 0, dipolar JJi] and supercooled 
liquids pol | , chemically reactive binary mixtures j2ll | . 

The dynamical equations will be solved using a nu- 
merical scheme based on a Lattice Boltzmann Method 
(LBM) for the Navier-Stokes equation and finite differ- 
ence methods for the convection-diffusion equation. LBM 
solves numerically a minimal Boltzmann equation where 
the fluid can only move along the links of a regular lattice 
with dynamics consisting of a free-streaming and a colli- 
sion step 22] . LBM has been largely applied to the s tudy 
of binary mixtures and complex flows [isL I23L l2^ I25I l26j| . 
Rheological behavior of liquid crystals under shear has 
been studied by LBM in Ref. [23. The mixed method 
used in this paper |2^ |2^ js^l allows a reduction in mem- 
ory which is convenient in large scale simulations. More- 
over, spurious terms appearing in Eq. (jSJ 24] are avoided 
and the numerical efficiency is increased. 

In our simulations we start from a disordered config- 
uration and consider the evolution of the system with 
parameters corresponding to the lamellar phase. This 
correspond to a sudden quench; shear is applied for all 
the evolution after the quench. Lamellae are expected 
to align with the flow with order propagating from the 
moving walls ^31j. One could think that after a suffi- 
cient long time lamellae will be ordered uniformly in the 
whole system. However, we will see that the system can 
behave differently. In some range of parameters we will 
find that a SIS region develops in the middle of the sys- 
tem consisting of small droplets and pieces of bent and 
rolled lamellae. This region coexists with lamellar re- 
gions close to the walls and a banded flow with different 
shear rates is observed. These results will be discussed 
in relation with the traditional scenario of shear banding 



the equilibrium model and the dynamical equations will 
be first introduced. Then the various components of the 
stress and the structure factor will be defined. Finally, 
a short review of the numerical scheme we use will be 
given. Section III contains our results for the evolution 
of the system in a specific case where a SIS region with 
banded flow is observed. The effects of changing viscosity 
and shear rate on the properties of the SIS region will be 
shown in Section IV. In Section V we analyze the evolu- 
tion of global order in the system in the spirit of what 
usually done in phase separation studies, considering the 
behavior of the first momenta of the structure factor. A 
discussion with conclusions will be given in Section VI. 



II. THE MODEL AND THE METHOD 

Our simulations are based on a mixed numerical ap- 
proach which combines the lattice Boltzmann method 
with a finite difference scheme H^H^. In this scheme the 
equilibrium properties of the system can be controlled by 
introducing a free energy. 



A. Equilibrium properties and dynamical equations 

The equilibrium phase is described by a coarse grained 
free energy that, for the specific problem we study here, 
is the following: 
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where n is the total density of the system and is a 
scalar order parameter representing the concentration 
difference between the two components of the mixture. 
The term in n gives rise to a positive background pres- 
sure and does not affect the phase behavior. The terms in 
If correspond to the Brazovskii free energy ^| . We take 
6, 0? > to ensure stability. For a > the fluid is disor- 
dered; for a < and k > two homogeneous phases with 
if = -iiyj —a/h coexist. A negative k favors the presence 
of interfaces and a transition into a lamellar phase can 
occur. In single mode approximation, assuming a proflle 
like A sin k^x for the direction transversal to the lamel- 
lae, one finds the transition (|a| — h) aX a —l.llK^/d 
where = ^-K/2d and =A{1 + K^/Mb)/3 [s^. 

The evolution of the system is described by a set of two 
coupled partial differential equations: The Navier-Stokes 
and the convection-diffusion equations. The fiuid local 
velocity u obeys, by assuming incomprcssibility (V • u = 
0) , the Navier-Stokes equation which reads as 



(2) 



where v is the kinematic viscosity. P^|^ is the thermody- 



The paper is organized as follows. In the next Section namic pressure tensor which can be calculated from the 
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free energy functional ^ as 

( 6F 6F 1 
P'^p^{n— + ^—-f{n,^)Yc.p + D^p{^) (3) 

where /(fi, is the free-energy density and a symmetric 
tensor Daf3{yy) has to be added to ensure that the con- 
dition of mechanical equilibrium daP^p = is satisfied 
[33|. The complete expression of the pressure tensor is 



p: 



af} 



(4) 



with 



1 



(5) 



and 



Moreover, we shear the system by moving the upper and 
lower walls with velocities 



L. 

±73!, 



(7) 



respectively, where 7 is the shear rate, L is the width of 
the system and i is a unit vector along the a;-axis, which 
is usually denoted as the flow direction. The presence of 
the moving walls greatly affects, as we will see, the fluid 
velocity u and the behavior of the order parameter Lp. 

The evolution of ip is described by the convection- 
diffusion equation: 



where 



6F_ 

5ip 



(8) 



(9) 



is the chemical potential [33|, F is a mobility coefficient, 
and the laplacian in the r.h.s. of © guarantees the con- 
servation of (fi. 



B. The shear stress 

The presence of shear strongly influences the morphol- 
ogy of the system and its rhcological properties. In par- 
ticular, we will consider the effects on shear stress. 

The total stress is (see the r.h.s. of Eq. 



Sal3 = -P^'p + nv {dfjUa + daUp) 



(10) 



The total shear stress, which is related to the off-diagonal 
part of the total stress 1(1^ , is 
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+ nv {dyUx + dxUy) 



(11) 



Sxy is the sum of the time reversible shear stress S^^™ = 
~Pxy'^"^ and of a dissipative hydrodynamic contribution 
5^^*' = nv{dyUx 4- dxUy). The expression (fTT|l depends 
on local coordinates and may vary from point to point 
especially when the system is not homogeneous. 

A quantity of relevant experimental interest is the in- 
tegral 

j dvSxy = S''^'-'"' + S^y'^'' = j drSfy''"''+ J drS^f"- 

(12) 

By using Eq. 10 it can be shown that S'^^'^"^ can be also 
calculated in the reciprocal space as 



gchera ^ , ^^g. 



dk 



-^—^kxkyiK + 2dk')C{k, t), 

(13) 

where C(k, t) is the structure factor 



C(k,t) = (^(k,i)^(-k,i)), 



(14) 



(y3(k) is the Fourier transform of the order parameter ip 
and (•) is the average over different histories. 



C. The numerical scheme 

The Eqs. (|2I8|I are numerically solved in 2D by using a 
mixed approach. We use the lattice Boltzmann method 
for Eq. and a finite difference scheme for Eq. ((HJ. 
Such an approach has been already adopted in the case 
of thermal lattice Boltzmann models for a single fluid and 
for multiphase flows ^231. In that case it is the tempera- 
ture equation to be solved by finite differences. By using 
this approach we are able to avoid the spurious terms 
in the convection-diffusion equation (jHJ which come into 
play when standard LBM for binary mixtures is used 0| , 
though LBM may realize boundary conditions easily and 
give better numerical stability. The present method al- 
lows to reduce required memory. 

A set of distribution functions /i(r, i) is defined on 
each lattice site r at each time t. Each function is 
associated to a lattice speed vector with e^/c = 
(±1,0),(0,±1),(±1,±1),(0,0), where c = Ax/ At, At 
is the time step, and Ax is the lattice constant. 

They evolve according to a single relaxation-time 
Boltzmann equation p3 . Issf : 

Mr + e,At, t + At)~ /,(r, t) - --[/,(r, t) - f^\r, t)] 

(15) 

where r is a relaxation parameter and f^''(r,t) are lo- 
cal equilibrium distribution functions. They are related 
to the total density n and to the fluid momentum nu 
through 



(16) 
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These quantities are locally conserved in any collision 
process and, therefore, we require that the local equilib- 
rium distribution functions fulfil the equations 

i i 

Following Ref . , the higher moments of the local equi- 
librium distribution functions are defined so that the 
Navier-Stokes equation can be obtained and the equi- 
librium thermodynamic properties of the system can be 
controlled: 

E f!''e^ae^p = C^P^^ + nUa^Up. (18) 
i 

The local equilibrium distribution functions can be ex- 
pressed as an expansion at the second order in the veloc- 
ity u 2^ ■ The expression of the coefficients of the equi- 
librium distribution functions can be found in Ref. |36| . 

The above described lattice Boltzmann scheme simu- 
lates at second order in At the continuity and the incom- 
pressible Navier-Stokes equations ^ with the kinematic 
viscosity v given by 

^ = At^(r-i). (19) 

It appears that the relaxation parameter t can be used 
to tune independently the viscosity. 

Equation (|SJ is numerically integrated by using an ex- 
plicit Euler scheme on a square lattice with spacing Ax, 
the same as for LBM. The spatial derivatives are approx- 
imated by discrete expressions which are second order in 
Aa;. The time step is At' — At/m with m = 5. This 
choice was motivated by the observation of a better nu- 
merical stability of the code. 

To enforce the flow Q we assume periodic bound- 
ary conditions (BC) along the flow direction and we 
place walls at the upper and lower rows of the lattice 
moving them at a constant velocity along the x direc- 
tion avoiding slip velocity [2^. The velocity u obtained 
from Eq. H16|) goes inside the convection-diffusion equa- 
tion ©. For the order parameter ip we adopt Lees- 
Edwards BC along the y direction |37| . This means that 
(f{x,-L/2) = ip{x + {vu - vi)At',L/2). The algorithm 
implementing the previous numerical scheme has been 
described in Ref. |28| . 

All the simulations in the following have been run by 
using the parameters —a = b = 2 x 10~^, k = —6 x 10~^, 
d = 7.6 X 10^"*, and F = 25. The system size was L = 
1024 and space and time steps were set to Ax = 1 and 
At = 0.2, respectively. At the beginning of each run the 
values of (p are randomly taken in the range [—0.1,0.1], 
the distribution functions fi are set so that n — 1, and 
the velocities u arc computed from Eq. (|16l) . We verified 
that by vertically shifting the reference frame, the fluid 
velocity profile is accordingly shifted. Therefore we used 
the choice of placing the walls symmetrically located at 
y = ±L/2 to have fiuid almost at rest in the middle of 
the system. 



III. KINETICS OF SIS FORMATION 

In this Section we will show the evolution of the system 
described by Eqs. H2I8|I for a typical case where shear 
induced structures appear. Shear rate and viscosity are 
fixed to 7 = 5 X 10~^ and v — 4.1667. All the quantities 
here and in the following are measured in units of the 
space step Ax and the time step At. For these values, 
the relaxation time of the linear shear velocity profile in 
a simple fluid would be of the order of t = 6000 [3^ . 

Figure ^ shows the configurations of (p at successive 
times. At t = 1000 after the quench, the order parameter 
has locally reached one of the minima of the polynomial 
part of the free-energy, represented by black and white 
in the figure. Lamellae are ordered only on small scales 
and in most of the system the effects of shear flow are not 
observable. Ordered structures appear only very close to 
the walls. At the next time t = 19000, lamellar order 
has developed into the system, but the the middle region 
is still almost isotropic apparently scarcely influenced by 
the flow. After this time the region of lamellae aligned 
with the flow does not increase the extension towards the 
central part of the system where morphology evolves in 
a different way. At i = 42000 the middle region mostly 
consists of lamellae oriented at about 45° with respect 
to the direction of the flow. Few domains are broken 
into droplets and small pieces of lamellae. The interface 
between the central and the two external lamellar regions 
is not very sharp and shows some undulations. 

The further evolution of the middle region can be seen 
a,t t = 60200. We will call this region as a SIS region or 
SIS phase. Many ruptures have occurred in the central 
network with the consequent formation of more droplets 
and worm-like domains. Not relevant changes can be ob- 
served at t = 81800 and t = 101800. Droplets, once 
formed, are quite stable. We checked the existence of the 
SIS region until time t = 200000. We repeated this nu- 
merical experiment starting from 5 different initial con- 
figurations and obtaining very similar results for the dif- 
ferent histories. 

The morphological evolution of the system has to be 
examined also in relation with the behavior of the veloc- 
ity profile. In Fig. |2 the x-averaged horizontal velocities 
Ux are plotted at the same times of Fig. ^ ^ functions 
of y. The variance corresponding to each average is also 
shown; it is generally quite small indicating an uniform 
behavior of the system in the flow direction. At the be- 
ginning (t — 1000) the shear proflle is different from zero 
only in a region of about 150 lattice sizes close to the 
walls; this explains the isotropy of the configuration ob- 
served almost everywhere at this time in Fig. ^ At later 
times the M^j-profile becomes characterized by two slopes 
found in correspondence of the regions with lamellar and 
SIS phases. The profile remains almost stationary from 
t ~ 40000; the local shear rates in the lamellar and SIS re- 
gions are respectively higher and lower than the imposed 
value. 

Mechanical properties are described by the behavior of 
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t = 1000 t = 19000 T = 4200C 




i = 6D200 i = SISOO t = lOlSOO 




FIG. 1: Configurations of tlie order parameter (p on tlie wliole 1024 x 1024 lattice for the case with 7 = 5 x 10 ^ and u — 4.1667. 



the two stress S^y'^™ and S!^!^'^^ . The structures present 
in the system locally determine the value of the chem- 
ical part of the stress. is larger in presence of 
defects or configurations which are not minima of the 
free-energy while it is close to zero for a well ordered 
lamellar configuration. This can be seen in the central 
part of Fig. 13 where we compare the chemical stress at 
different distances from the walls. In the middle of the 
system (y = 0) the x-averaged stress reaches a maximum 
in correspondence of the formation of the SIS and then 
remains different from zero. Fluctuations correspond to 
the evolution of structures present in the SIS phase. The 
behavior of the chemical stress aX y = L/8, just inside 
the SIS region, is similar. On the other hand, in the 
region with lamellar order, at y = L/A and y = 3L/8, 
after a maximum at initial times, the stress relaxes to 
a very small value. This maximum occurs before than 
lamellae become aligned with the flow. This behavior is 
analogous to that observed in phase separation of binary 
mixtures under shear where the excess viscosity reaches a 
maximum before than domains orientate with interfaces 
in the direction of the flow (39|. 



The hydrodynamical stress, also shown in Fig. O has 
an opposite behavior: It is larger where shear rate is 
larger (y = i/4 and y — 3L/8). Closer to the walls 
{y = 3L/8) it relaxes to a constant value after an initial 
maximum corresponding to the fact that the shear rate 
close to the walls is higher at initial times. The bottom 
part of Fig. |31 shows that the system is in mechanical 
equilibrium from t ^ 20000 onwards when the total stress 
has become the same at different distances from the walls. 

Finally, in Fig. 21 the global patterns of the velocity 
field are shown for the same times of Fig. ^ We consider 
a central region of the system which includes a portion 
of the SIS phase and the interface with the lamellar re- 
gion. At t = 1000 there is no evidence of the imposed 
fiow in this part of the system; local flows, as usually in 
phase segregation, originate from the defects. Effects of 
the imposed fiow becomes evident a.t t = 19000. Moving 
inside the SIS phase, the magnitude of the velocity vector 
becomes smaller with horizontal and vertical components 
comparable. Vortex structures, not present in the lamel- 
lar phase, can be also observed. In later pictures one 
observes that the SIS-lamellar interface structure evolves 
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FIG. 2: Plots of the a;-averaged horizontal velocities across 
the system at different times for the case with 7 = 5 x 10~* 
and V = 4.f667. 



with time together with the local flow patterns. Veloc- 
ity in the SIS phase always remains small and a jump in 
the velocity magnitude can be seen moving through the 
SIS-lamellar interface at all the late times considered. 
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IV. PROPERTIES OF SIS PHASE 

The behavior described in the previous Section is 
found, at a given viscosity, for shear rates under a cer- 
tain threshold. Velocity, stress and other properties, in 
cases when a SIS region is observed, have been checked 
to remain stable for long times and a linear velocity pro- 
file is never reached. Figure |31 shows the behavior of the 
chemical part of the stress for two values of shear rate 
smaller and larger than the threshold that, for the case 
V = 4.1667, is 2.5 x 10"'' < 7c < 5 x lO^^. 

At 7 = 10"'* the chemical stress remains almost sta- 
tionary after the formation of the SIS region and the mor- 
phological evolution is similar to that shown in Figs. 11121 
On the other hand, at 7 = 10~^, after a maximum, the 
shear stress decreases to small values while the horizon- 
tal velocity tends to a profile with constant slope. In this 
case the SIS region does not form in the middle of the 
system and ruptures do not occur at all in the lamellar 
network. Late-time configurations for the two cases are 
shown in Fig. El For 7 = 10"'^ one sees a well-ordered 



FIG. 3: Plots of the different components of stress, summed 
over the a;-axis, as functions of time at different positions 
across the system for the case with 7 = 5 x 10~^ and 
u = 4.1667. 



lamellar state with only few local defects. 

For the case v = 4.1667 we run various simulations 
with different values of 7 < 7c also for comparing our 
results for the SIS pha se with the usual scenario of shear 
banding Ullilllli^. In this scenario bands of phases 
with different structure and flow properties can coexist 
only at a critical value of the stress S* . In experiments 
with imposed shear rate this leads to a stress plateau on 
the flow curve (stress versus imposed shear rate) at S* . 
On the plateau, the shear rates of the different bands 
do not change by varying the imposed flow which only 
determines the relative spatial extension of the bands. 

In Fig. [7| we show the shear rate profiles for different 
values of 7 < 7c and for 7 = 5 x 10~^ at times, for 
7 < 7c, when the SIS phase is well formed. A central 
region with smaller shear rate can be clearly identified in 
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FIG. 4: Configurations of tiie order parameter <^ on a central 
portion of size 128 x 128 of the wiiole lattice for the case with 
7 = 5 X 10~^ and v — 4.1667. Arrows represent fluid velocity 
vectors u and are plotted every four lattice sites. 



all cases, except that for7=:5xlO^^. The width of this 
region, the one with SIS morphology, decreases when the 
imposed sear rate is increased, as usually in systems with 
shear banding. Close to the walls the shear rate is larger 
than the imposed value with an almost constant value in 
correspondence of the region with well aligned lamellae. 
However, differently from other systems with shear band- 
ing, the values of shear rates in the lamellar and in the 
SIS regions change with the imposed flow. The constraint 
that the integral of the horizontal velocity in the verti- 
cal direction gives the value imposed on the boundaries 
is always verified. An interface of finite width between 
the lamellar and the SIS region can be seen in all cases, 
confirming the relevance of the arguments discussed in 
Ref. m. 

The behavior of the hydrodynamical, chemical and to- 
tal stress is shown in Fig.|Slas a function of the imposed 
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FIG. 5: Plots of the integral chemical stress (Eq. as a 

function of time for two cases with different shear rates. 



shear rate. The total stress does not exhibit a flat regime 
in the interval where the SIS phase exists and changes in 
the behavior are not observed when 7 becomes greater 
than 7c. As expected, the chemical stress decreases with 
the reduction of the width of the SIS region. 

Finally, we examine the effects of viscosity showing re- 
sults at fixed shear rate for times after the SIS formation. 
The main result is that at higher viscosities the width of 
the SIS region becomes narrower, as shown in Fig. |51 
This behavior is confirmed from Fig. IIUI where shear rate 
profiles arc plotted. 

We have tried to interpret these results considering 
possible relations between cases with a SIS region and 
characteristic time scales of the imposed flow. The shear 
rate naturally introduces a typical time — 7"^ when 
shear effects are expected to become relevant 0. On 
the other hand, the propagation of the horizontal velocity 
from the moving walls requires a relaxation time which, 
in simple fluids, is inversely proportional to the viscosity 
[ssj l . This time is also relevant for the effectiveness of the 
imposed flow in the middle of the system. We did not 
find a quantitative explanation of the phenomena before 
described in terms of these time scales, probably, also 
because the relaxation time of a simple fluid is not ap- 
propriate for the system we are studying. However, our 
results suggest the following general qualitative observa- 
tion. If shear effects arrive early enough in the middle 
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(a) 





FIG. 6: Late time configurations of the order parameter Lp on the whole 1024 x 1024 lattice with v = 4.1667 in the cases with 
7 = 10"^ at time t = 26400 (a) and 7 = 10"" at time t = 73900 (b). 
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FIG. 8: Plots of the different components of stress summed 
over the whole system at long times as functions of the im- 
posed shear rate 7 for the case with v — 4.1667. In the inset 
the chemical stress is plotted on a magnified scale for a better 
view. 



200 400 600 



-600 -400 -200 



200 400 600 



FIG. 7: Plots of the a;-averaged measured shear rate across the 
system for different imposed shear rates 7 with f = 4.1667. 
Values are normalized to the imposed shear rate 7. 



of the system, due to a large shear rate or to a high vis- 
cosity, the flow is able to penetrate in this region and 
lamellae align with it. Otherwise, at lower shear rates or 
lower viscosity, a SIS region can be found. We will return 
on this point in the discussion of the last Section. 



V. DOMAIN GROWTH PROPERTIES 

In cases without imposed flows, the kinetics of forma- 
tion of ordered phases after a quench is a process which 
is well-known for simple systems like binary mixtures, 
and still under investigations for lamellar systems. In 
binary mixtures dynamical scaling occurs with domains 
growing isotropically with power-law behavior ^41] . In 
lamellar systems the formation of extended defects or 
grain boundaries between domains of differently oriented 
lamellae characterizes the late time evolution; other de- 
fects like dislocations and disclinations are also present 
[sol li^ li^ . Typical lengths can be defined and their 
evolution follow power-law or slower asymptotic growth 
depending on the steepness of the potential energy and on 
the depth of the quench. In cases with power-law behav- 
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FIG. 10: Plots of the a;-averaged measured shear rate across 
the system for different viscosities v and 7 = 10"''. 




ior the length corresponding to the inverse of the struc- 
ture factor at half height grows with an exponent ranging 
in the interval 0.2 < z < 0.33 [H El H gl El E El ■ 
In our case with shear the presence of banded struc- 
tures makes problematic the description of ordering in 
terms of a single length, as we have seen in the previous 
sections. However, it could still be interesting to consider 
the behavior of averaged lengths intended as macroscopic 
indicators of the degree of order in the system. We define 



/ dkC(k,t) 
/ dk|fc„|C(k,i) 



a ■ 



(20) 



and show in Fig. ^2 the evolution of Rx,Ry for differ- 
ent viscosities and shear rates. In all cases Ry relaxes to 
the equilibrium value, while R^ grows with exponents 
varying with the viscosity and the shear rate as reported 
in the figure. In particular, keeping the viscosity v fixed, 
we found that the value of Zx increases by increasing the 



FIG. 11: Plots of the domains sizes Rx and Ry as functions 
of time for different values of viscosities and imposed shear 
rates. See the single frames for parameters details and slopes 
of straight lines. 



shear rate 7. In the case discussed in Section III with 
7 = 5 X 10"^ it results that z^ = 0.38 for times t < 40000 
corresponding to SIS formation; the later growth is char- 
acterized by the elimination of local defects in the lamel- 
lar region. Finite size effects may affect the behavior of 
Rx at very late times. The case not shown in Fig. ITTI 
with 7 = 10~^ and the same viscosity v — 4.1667 is simi- 
lar. The time interval corresponding to the SIS formation 
{t < 20000) is characterized by a growth compatible with 
exponent z^ — 0.52 while the growth rate changes later. 
By increasing further the shear rate (7 = 10"'^), banded 
configurations do not form (see Figs. I5I6|I and the whole 
system reaches a lamellar order with exponent z^ ^ 1.2. 
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The behavior at 7 = 10~^ and v = 7.5 is quahtatively 
similar to those aheady discussed; the difference is in the 
value of the growth rate which is now z^. = 0.6 during the 
SIS formation. We observe that this value is greater than 
in the case with the lower viscosity v — 4.1667. Finally, 
for the highest value of viscosity shown in Fig. 1111 a power 
law appears not appropriate for describing the behavior 
of Rx ■ We can conclude that the growth of lamellar order 
as measured by Rx is faster for parameters correspond- 
ing to a smaller extension of the SIS region and a larger 
lamellar phase. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have shown results from two- 
dimensional simulation results for a fluid mixture 
quenched from a disordered configuration into a state 
with lamellar order. The fluid, described by Navier- 
Stokes and convection-diffusion equations, is subject dur- 
ing all the evolution to the action of a shear flow imposed 
by the walls of the system. The ordering process of a 
lamellar system, in absence of flow, would be character- 
ized by the presence of local and extended defects that 
make slow, sometimes freezing, the evolution of the sys- 
tem. When shear is applied, defects tend to be eliminated 
and lamellae would be expected to align in a preferred di- 
rection forming well-ordered macroscopic domains. 

Our results actually show that different evolutions with 
more complex morphologies are also possible under shear. 
For small enough shear rates and viscosities we found 
that the flow stabilizes itself with a horizontal velocity 
profile characterized by two different slopes. Close to 
the walls the shear rate is higher than the imposed one 
and well-ordered lamellae can be observed. In the central 
part of the system, that is reached later by the flow, the 
shear rate is smaller than the imposed one and the mor- 
phology of domains is characterized by the presence of 
small droplets and pieces of bent or rolled lamellae never 
aligned with the flow. We referred to this region as to a 
SIS phase. 

Shear banding phenomena occur in many complex flu- 
ids. Similarly to what generally found, in our case the 
width of the SIS phase decreases by increasing the value 
of the imposed shear rate. In other important aspects, 
however, our results differ from the usual picture of shear 
banding. In the range of shear rates with coexisting 
lamellar and SIS regions, we do not observe a plateau 
of the total stress at varying shear intensity. Moreover, 
in this range, the values of local shear rates corresponding 
to the two phases depend on the imposed flow. 

At the moment, we are not aware of experiments with 
the same scenario as that described in this paper. How- 



ever, even if our model can be appropriately used for 
studying copolymer systems in the weak segregation limit 
|16| | , the comparison with experiments is limited from the 
fact that our simulations are two-dimensional while in 
three-dimensional systems more complex geometries can 
occur [iH I49I I , with possible different rhcological behav- 
iors and flow patterns. Actually, the purpose of this work 
was more generic. We wanted to analyze general features 
of the formation of banded flows in systems with lamel- 
lar order considering the full dynamical problem for the 
velocity and the concentration fields. 

Our results of Section III and IV show that when the 
applied flow is weak enough or propagates from the walls 
sufhciently slowly (at lower viscosities), it is not able to 
penetrate with the same shear rate in all the system. In 
this case the central region, with its intertwined tangled 
lamellae, opposes to the presence of the flow and tends 
to keep its morphology. The region tends to behave as 
an almost frozen network; the shear mainly acts stressing 
this network and causing ruptures with the production of 
droplets and small pieces of lamellae. Even if the shear 
rate in the SIS region is quite small but not zero, the be- 
havior of our system resembles that of yield stress fluids 
where no flowing steady states exist for stresses under 
a certain threshold For example, in soft glass sys- 

tems, coexistence between not flowing "pasty" states and 
sheared fluid regions has been observed In rela- 

tion with this, it is interesting to observe that an equi- 
librium glass transition has been found also in systems 
with lamellar order j^]. We think that in our case the 
competition between the intrinsic slow dynamics [8(ij | and 
the acceleration induced by the external flow is respon- 
sible for the peculiar shear banding phenomena we have 
shown. 

Finally, we observe that our numerical methods have 
been proven to be quite convenient for simulating fluids 
with complex order under driving forces. These methods 
can be implemented also with different geometries and 
in the more realistic three-dimensional case. We hope 
that our two-dimensional results can be useful for the 
comprehension of shear banding and stimulating for new 
experiments. 
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